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ABSTRACT 

We propose that the internal energy of the GRB blast waves, thought to be 
stored in the form of relativistic protons co-moving with the blast wave, is con- 
verted explosively (i.e. on light crossing time scales) into relativistic electrons 
of the same Lorentz factor, which are responsible for the production of observed 
prompt 7-ray emission of the burst. This conversion is the result of the combined 
effects of the reflection of photons produced within the flow by upstream located 
matter, their re-interception by the blast wave and their eventual conversion into 
e"'"e~— pairs in interactions with the relativistic protons of the blast wave (via 
the p7 — > e+e" reaction). This entire procedure is contingent on two conditions 
on the relativistic protons: a kinematic one imposed by the threshold of the 
P7 e^e~ reaction and a dynamic one related to the column density of the post 
shock matter to the same process. This latter condition is in essence identical 
to that of the criticality of a nuclear pile, hence the terminology. It is argued 
that the properties of relativistic blast waves operating under these conditions 
are consistent with GRB phenomenology, including the recently found correla- 
tion between quiescence periods and subsequent flare fiuence. Furthermore, it 
is shown that, when operating near threshold, the resulting GRB spectrum pro- 
duces its peak luminosity at an energy (in the lab frame) E ~ nieC^, thereby 
providing an answer to this outstanding question of GRBs. 

Subject headings: accretion, accretion disks — radiative transfer — line: forma- 
tion — X-rays: general 
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1. Introduction 

The longstanding issue of the distance and absolute luminosity of GRBs has been settled 
in the past decade through the observational evidence collected by BATSE (Meegan et al. 
1992) and BeppoSAX (Costa et al. 1997), while the theoretical work of Meszaros & Rees 
(1992) and Rees & Meszaros (1992) provided the broader physical framework into which 
these events seem to generally fit. This framework associates GRBs with radiation emitted 
by relativistic blast waves (hereafter RBWs), produced by an unspecified todate agent, 
presumably associated with the formation of a neutron star or black hole. While the source 
of the energy associated with GRBs has remained uncertain, there remains little doubt about 
the presence of the RBWs, which power also the later time emissions at X-ray, optical and 
radio frequencies, known collectively as GRB afterglows (see review of Piran 1999). 

With the discovery of GRB afterglows, much of the theoretical activity has since shifted 
to the study of the physics of these later time emissions. Nonetheless, a number of issues 
associated with the prompt 7-ray emission, besides the nature of their energy source, still 
remain open. Chief among them are the conversion of the RBW kinetic energy into radiation 
and the fact that the frequency at which the GRB luminosity peaks, Ep, is narrowly dis- 
tributed around a value intriguingly close to the electron rest mass energy. The purpose of 
the present note is to describe a process that provides a "natural" account of these generic, 
puzzling GRB features. 

Following the work of Shemi & Piran (1990) it has been generally accepted that a certain 
amount of baryons must be carried off with the blast waves responsible for the GRBs. This 
baryon contamination has even been deemed necessary, else the entire blast wave internal 
energy would be converted into radiation on very short time scales, leading to events of very 
different temporal and spectral appearance (e.g. Paczyhski 1986) than observed in GRBs 
. While the low radiative efficiency of baryons is essential for the GRB energy transport 
to the requisite distances ( ^ 10^^ cm), it becomes problematic when demanded that their 
internal energy be radiated away on the short time scales associated with the GRB prompt 
7-ray emission. Generally, this issue is sidestepped by appealing to an unknown process 
which transfers the proton energy into electrons (Dermer, Bottcher & Chiang 1999), whose 
radiative evolution could then be accurately computed. 

The narrow range of the GRB uFi, spectral peak energy, Ep, is another well docu- 
mented systematic feature of these events, a result of the extensive data base accumulated 
by BATSE. The compilation of Malozzi et al. (1995) shows clearly a preference for an 
energy Ep ~ 200 keV at which the uF,^ GRB spectra exhibit a maximum. In fact, it is 
precisely this maximum in the spectral energy distribution that qualifies GRBs as such. 
Furthermore, when corrected for the redshift (-^grb ~ 1), Ep shifts close to the electron rest 
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mass. While a compelling explanation of this fact is presently lacking, several accounts have 
occasionally been proposed. For example, Brainerd (1994) argues that this is the result of 
down-Comptonization of a power law photon distribution that extends to E ^ 1 McV by 
cold matter with Thompson depth ~ 10, an explanation possibly in conflict with the 
timing properties of GRBs (see e.g. Kazanas, Titarchuk and Hua 1997). The association of 
the GRB emission with relativistically boosted synchrotron radiation from RBWs has made 
this particular issue far more acute, as the energy of the latter should scale like Ep oc F^. 
Therefore, even very small variations in the values of F would lead to a very broad range in 
the values of Ep. Dermer et al. (1999) proposed that the observed distribution is the result 
of the time evolution of a blast wave with a specific baryon loading, which when convolved 
with the triggering criteria of existing detectors favors the detection of fireballs with Ep in 
the observed range. On the other hand, on the basis of analysis of SMM data, Harris & 
Share (1999) have argued that there is no apparent excess of GRBs with Ep ^ 1 MeV, thus 
leaving this issue open. 

The present paper is structured as follows: In §2 we outline the fundamental notion 
behind our proposal for converting the RBW proton energy into radiation, we derive the 
associated necessary conditions and discuss its relation to GRB phenomenology. In §3 we 
produce model spectra based on this proposal and indicate their relation to the particular 
value of Ep observed. Finally, in §4 the results are discussed and certain conclusions are 
drawn. 



2. The "Supercritical Pile" 

The process described herein has been discussed in the past by Kazanas & Mastichiadis 
(1999; hereafter KM99) in the context of AGN, where arguments have been put forward in 
favor of a hadronic origin of the relativistic electrons in blazars. This process is effectively 
the relativistic plasma instabihty proposed by Kirk & Mastichiadis (1992; hereafter KM92) 
coupled to the increase in the photon energy associated with relativistically moving "mirrors" 
(Ghisellini Sz Madau 1997). While the mathematical formulation of the instabihty is given 
in detail in the above references, we provide below a qualitative re-derivation which focuses 
on and elucidates the underlying physics. 
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2.1. 



A Static Plasma 



Consider a spherical volume of size R, containing a relativistic proton plasma of dif- 
ferential spectrum 72^(7) = no7~'^ (7 being the proton Lorentz factor), along with an (in- 
finitesimal) number of photons of energy e (in units of mgC^); these photons can produce 
pairs via the p7 — > e~e+ reaction, provided that the proton population extends to Lorentz 
factors 7 > 7c such that 7c e ~ 2. In the presence of a magnetic field B, the pairs (of Lorentz 
factor also equal to 7c) produce synchrotron photons of energy eg = 67^ where b — B/Bcr is 
the magnetic field in units of the critical one Bcr — n?^(? j (eh) ~ 4.4 10^^ G. 

For the reaction network to be self-contained the energies of the seed and synchrotron 
photons should be equal, yielding the kinematic threshold of the process i.e. 7ces = 7^ 6 ~ 2. 
For the process to be also self-sustained, at least one of the synchrotron photons must produce 
a pair before exiting the volume of the source. Hence, the optical depth of the source to 
the p7 — s> e~e~^ reaction should be greater than l/A/", where A/" ^ 70/67^ = l/&7c is the 
total number of synchrotron photons produced by an electron of energy 7c. This condition 
then reads Tp^ ~ Rnp{j)j = ap^ Rno%'^^^ ~ ^7c- Eliminating 7c using the kinematic 
threshold condition 7c 6 — 2, the critical column density expression reduces to 



which (within factors of order unity) is the condition derived in KM92. This condition is 
similar to that of the criticality of a nuclear pile, except that there one deals with neutrons 
rather than photons (this similarity carries also over in the case in the Comptonization of 
photons by hot electrons (Katz 1976)). It becomes apparent, hence, that the critical quantity 
here (as also in a nuclear pile) is the column density (rather than the mass, the term "critical 
mass" being simply a figure of speech). The pair - synchrotron photon - proton - pair network 
will be self-sustained if the column density is equal to the critical one. For larger values the 
number of pairs increases exponentially, eventually leading to a depletion of the available 
energy source on time scales ~ i?/c (a bomb!). 



KM99 extended the above analysis to the case that the relativistic proton containing 
plasma moves itself relativistically with Lorentz factor F. The criticality conditions change 
quantitatively if the radiation emitted by the plasma can be "refiected" by matter located 
along its direction of motion. Due to relativistic beaming, essentially all the photons pro- 
duced internally in the plasma will be focused in the forward direction, reflected and boosted 
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2.2. 



Plasma in Relativistic Motion 
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in energy, upon their re-interception by the moving plasma, by a factor F^. If 6 is the (nor- 
mahzed) value of the comoving magnetic field, the synchrotron photons of energy = 67^ 
will, upon their rc-intcrception, have energy = b^^F^, modifying the kinematic threshold 
condition to by^^ 

The change in the kinematic threshold affects also the dynamic one: The photons neces- 
sary for the production of pairs in the comoving frame arc now emitted by electrons of Lorentz 
factor only 7c/r. The number of such photons is now A/" — (7/r)/6(7/r)^ = l/6(7/r). De- 
manding again that the column of the plasma be greater than l/N', along with the new 
threshold relation, leads to the condition 

PI ^com "'Ocom ^ ^ V J 

namely the condition derived in KM99 (Acom, ''^ocom are the co- moving source size and den- 
sity). Relativistic motion therefore, eases significantly the "criticality" condition in the case 
of a relativistically moving plasma. 

The situation in the RBW of a GRB is somewhat different than that discussed just 
above. Though diffusive acceleration is very likely present in their associated shocks, the 
mere postshock isotropization of the flow creates, in the RBW frame, a relativistic population 
of protons with mean energy {E) ~ TmpC^. Therefore, to be most conservative, one could 
dispense with the requirement of an accelerated proton population and demand that only 
the protons of energy FmpC^ be present (these are certainly the most numerous). Therefore, 
upon setting 7^. ~ F, the kinematic threshold condition reduces to T^b ^ 2. The number 
of photons emitted by electrons at threshold then becomes H — 1/b and the "criticality" 
threshold reads 

^Vl ^com ^com ~ RtT- ^ b OT Up^ RflF ^2 (3) 

with the last relation incorporating both the kinematic and dynamic threshold (Acom '^com = 
Rn due to the Lorentz invariance of the column density. In this case, the quantity n R simply 
denotes the amount of ambient matter per unit area swept by the RBW). For the typical 
values of n and R used in association with GRBs, i.e. n = 1 Uq cm~^ and R = 10^^ Riq cm 
and considering that ~ 5 10~^^ cm^, the criticality condition yields F ^ 180 (uq Rig)~^^^, 
values well within the accepted parameter range. 

Following the above analysis, the correlation between quiescence and activity periods in 
GRBs found by Ramirez- Ruiz & Merloni (2000), finds a direct, qualitative interpretation: 
A RBW sweeps and 'piles-up' the ambient medium behind its forward shock, whose column 
density increases from below to above its critical value. Provided that the proper "mirror" 
is located upstream, the internal energy, or part of it depending on conditions, will be 
explosively converted into radiation. The energy released then would be proportional to the 
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amount of matter accumulated during quiescence, with the process repeating as the blast 
wave progresses and more matter is accumulated. 

3. The Spectra 

Generally, the hght curves as well as the GRB spectra are highly variable. In the present 
model, the issue of variability becomes even more complex than in more conventional models, 
given that the injection rate of relativistic pairs depends on proton collisions with photons 
emitted at prior times, following their reflection by the "mirror" . 

While the computation of the spectra within this model is an inherently time dependent 
problem, we present herein a simplified, steady-state treatment whose salient features, we 
believe, will be preserved in a more detailed calculation. 

We consider a RBW of energy E — 10^^ E^i erg in a uniform medium of density n—lriQ 
cm~^ and of a (normalized) opening half angle Q — (O/n). Because of the relativistic focusing 
of radiation, we need only consider a section of the blast wave of opening half angle 6 — l/F. 
The shocked electrons of the ambient medium (and pairs from the P7 — > e'^e~ process) 
produce synchrotron photons of energy e^. ~ bV'^. These, upon their scattering by the 
"mirror" and rc- interception by the RBW, are boosted to energy e = e^F^ = feF^ (in the 
RBW frame). These photons will then be scattered by (a) electrons of 7 ~ 1, originally 
contained in the RBW and/or cooled since the explosion, and (b) by the hot (7 ~ F), 
recently shocked ones to produce inverse Compton (IC) radiation at energies correspondingly 
€1 ~ 6 F^ and £2 ~ & F^ at the RBW frame. At the lab frame, the energies of these three 
components, i.e. e^, ei, €2 will be higher by roughly a factor F, i.e. they will be respectively at 
energies b F^, b F^ and b T^. Assuming that the process operates near its kinematic threshold, 
6F^ ~ 2, at the lab frame these components will be at energies ~ F~^, ei ~ 2 ~ 1 MeV 
and 62 — F^ ~ 10 GeV (F/100)^. This model therefore, produces "naturally" a component 
in the uF^, spectral distribution which peaks in the correct energy range. It also predicts the 
existence of another component at an energy F^ higher; such high energy emission has been 
observed from several GRBs (Dingus 1995). 

The relative importance of these components depends on the specifics of a given source. 
For the two IC components this depends on the ratio of the scattering depths of the cool 
rT(7 ~ 1) and hot ttI'J — F) electrons. At a minimum, this ratio is equal to the ratio 
of their escape to cooling rates T^-cooi/'^esc — ^esc/^cooi- li Ub = f b{E / At^ Q"^) is the 
magnetic energy density (/s denotes departures from equipartition; ~ 1 in GRBs), then 
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Fig. 1. — The inferred isotropic luminosity of the three components discussed in the text for 
a burst with AQ/47r = 0.005, r= 240 and kinetic luminosity Lt ~ 2 x lO^'^ erg/s. The dot- 
dashed line is the synchrotron component, while the solid and dashed lines are respectively 
the Compton scattered synchrotron by the cold and hot electrons. The dotted horizontal 
line denotes the equivalent isotropic kinetic luminosity. For a duration of ~ 20 s. this will 
result in a fluence of 10^^ ergs. 



tesc — Acom/c ^ R/T c and ^cooi — rUeC^/ {Ub (7t cF) so that 

Tr(7~l) fsfcEsi 



(4) 
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where /c(> 1) is a multiplicative factor indicating the contribution of cold electrons inherent 
in the RBW to the scattering depth. For G ~ 1/F one obtains for the rate of photon 
scattering into the two IC components 

^ r^) ^ ^ 1) , (5) 

jBjci- -C/51 

indicating that their luminosities would be roughly equal for the fiducial values of the pa- 
rameters involved. 



The relative luminosity between the direct synchrotron and the cold-electron scattered 
synchrotron photons can be estimated as follows: The number of electrons swept-up to radius 
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by a section of transverse dimension dt ~ R/T of the RBW is Afe ~ nE? fm/^^, where 
/m(< mp/me) is a multiphcative factor denoting the number of pairs produced in the RBW. 
Given that each electron of energy F produces ~ l/6r photons, the total number of photons 
produced to radius R is 

. /- IT' fm R /„\ 

These photons will be received in the lab frame at energy ~ bV^ over a time interval 
~ R/r'^c to yield a photon flux 

• nfraR^c ^ f Rc\ 

M.c^-^-^c^nRa^r [— ) U (7) 

It is assumed that a fraction a of these photons will scatter at the "mirror" and randomize 
producing a "layer" of photons of width A ~ R/T"^ (in the lab frame; it is assumed that 
the "mirror" is thinner than -R/F^), density ~ Af^a/V — Afja/{R/rY A and Thompson 
depth tt,7 ^ {naTR/br)a fm (Tr,7 — na^ RT^ a fm if near threshold i.e. for hV ~ l/F^)- 

The number of electrons to radius R swept by the section of the RBW considered here is 
then Me — {nR^ /r^)(/c + fm)', hence the total number of photons scattered by the electrons 
upon traversing the "photon layer" of width A will be N'sc,j — ■A/'e 'Tt,^- These will be received 
at the lab frame over a time ~ A/cF^ ~ R/cV^ to yield a photon flux (at energy hV^ ~ 2) 

Nsc,^ ~ (ni?aT)'F6 afm{fc + fm) (8) 

suggesting that the ratio of the photon fluxes given by Eqs. (7), (8) will be proportional to 
nRa^r F^q;(/c + fm) ■ For reasonable values of the source parameters one does obtain similar 
photon fluxes at these two energy bands (6F^, 6F^) indicating that the luminosity at ~ 1 
MeV is roughly F^ times that at the optical band. 

In figure 1 we present a sample of a spectrum obtained using the arguments discussed 
above. Assuming the process to operate near threshold (with no accelerated proton compo- 
nent), leads to S— function like injection of pairs; this results in the vFi, oc i/^/^ spectrum 
of figure 1 (which at low energies reverts to the uFi, oc i/^/^ spectrum of thin synchrotron). 
Clearly more involved spectra will result if particle acceleration is also incorporated in these 
models; however, this will be the subject of future work. 

4. Discussion 

We have presented above a novel model for the prompt emission of GRBs. We believe 
that this model provides some of the missing physics between the RBW proposal, which 
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describes successfully the GRB energetics and time scales and the prompt emission of ra- 
diation (7-ray as well as optical-UV), by a well defined mechanism for tapping the kinetic 
energy stored in the RBW baryonic component. Furthermore, the same physics employed 
in effecting the conversion of baryon kinetic energy into radiation is instrumental in produc- 
ing a peak in the spectral energy distribution at £'p ^ 1 MeV, thus providing a "natural" 
account of this GRB feature, an issue that has actually gotten even more puzzling with the 
advent of the RBW model for GRB. In addition to this emission, the model implies the 
presence of 7-ray emission at higher energies, namely E ~ F^meC^, a fact supported by some 
observations already, but which will be explored in greater depth by SWIFT and GLAST, 
missions which will be able to test and put meaningful constraints on this specific model. 

One of the fundamental features of this model is the presence of the combination of 
kinematic and dynamic thresholds; this combination, along with the presence or absence of 
the necessary "mirror", make the model inherently time dependent. At the same time, the 
kinematic threshold provides (for the first time to our knowledge) a regulating mechanism 
that puts one of the peaks in the vF^ distribution very close to the observed value, despite 
the motion of the emitting plasma with Lorentz factors of several hundreds. Concerning the 
nature of the "mirror" we are willing to speculate here that, because the main "reflected" 
component consists of the prompt synchrotron photons which are emitted at optical frequen- 
cies, one could use atomic cross-sections (~ 10~^^ cm^) to estimate their reflected fraction 
a; values of n, it! typical to those associated with GRBs yield a ~ 0.01 — 1. This is a very 
rough estimate, because it ignores the ionization of the reflecting medium. A more detailed 
treatment including these effects is beyond the scope of the present work. 

We would like to thank Jay Norris, Rob Preece and Brad Shaffer for stimulating dis- 
cussions and much information on GRB phenomenology. 
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